Can annihilating Dark Matter be lighter than a few GeVs? 
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We estimate the gamma ray fluxes from the residual annihilations of Dark Matter particles having 
a mass rridm £ [MeV,0(GeV)] and compare them to observations. We find that particles lighter 
than 0(100 MeV) are excluded unless their cross section is S-wave suppressed. 



o 
o 

(N 

m 
> 
oo 
in 

00 

o 

O 

Oh 
I 

O 



Introduction 

The accurate measurement of galactic rotation curves, 
the CMB spectrum, the primordial abundances of light 
elements, together with our understanding of structure 
formation provide convincing evidence in favor of the 
existence of Dark Matter (DM) [3. While the MA- 
CHOs searches 2] indicate that an astrophysical solu- 
tion is rather unlikely, most efforts are now concen- 
trated on searches for Weakly Interacting Massive Par- 
ticles (WIMPs) 0. These particles would belong to the 
Cold Dark Matter scenario (CDM); they would annihi- 
late and suffer from negligible damping effects at a cos- 
mological scale. Considering fermions only and assum- 
ing Fermi interactions, it was concluded 4] that the relic 
density argument constrains the DM mass {m^m) to be 
greater than a few GeVs, as is quite naturally predicted 
within the framework of supersymmetry. Nevertheless 
searches for very massive particles {vridm ^ O(GeV)) re- 
main unsuccessful so there is still room for new sug- 
gestions. 

Alternative DM scenarios have been proposed in re- 
sponse to the discrepancy between observations and 
CDM numerical simulations on small scales (which no- 
tably predict cuspy haloes 0). The most "robust" 
one, the Warm Dark Matter scenario (WDM), involves 
non-annihilatinq particles and a very narrow range for 
the DM mass [7|, obtained by requiring that the free- 
streaming length be of the order of the smallest primor- 
dial scale one wants to be compatible with. For instance, 
a typical warmon mass is mdm ^ 0( keV) for ~ 100 kpc 
but this scenario of light and non-annihilating DM, 
although not excluded as yet, also gives cuspy haloes 0- 

In this letter we propose another model, somehow 
intermediate between CDM and WDM. We shall in- 
deed consider annihilating DM particles having a mass 
rridm G [MeV, 0( GeV)]. This has rarely been stud- 
ied, probably because of the Lee-Weinberg argument 
that we shall evade here by considering bosonic (in fact 
scalar) candidates (for mdm ^ GeV) . Some of these par- 
ticles could perhaps turn out to be Warm not because 
of their mass but because of their collisions with rela- 
tivistic species If ^.ot excluded by any cosmolog- 

ical/astrophysical arguments, they could compete with 
the collisionless WDM and CDM scenarios. On the 



other hand, they are likely to fail to give flat galactic 
cores at '--^ 1 kpc despite their quite large annihilation 
rate and a possible significant damping mass. Interest- 
ingly enough, they should escape present DM direct de- 
tection experiments (which so far are only sensitive to 
masses greater than ~ O(GeV)), as well as accelera- 
tor experiments, as briefly discussed in the next section. 
These particles would be compatible with the blackbody 
spectrum measurement and will not yield any *i?e pho- 
todissociation (for mdm > 26 MeV) provided their (s- 
wave) cross section satisfies the relation (mdm/ MeV) > 
5 [(crWr)arm/3 X 10"^'' cm^ s~^] (fldmh"^)^ (assuming self- 
conjugate DM particles, 7J = 100/ikm/ s/Mpc and using 
the D measurement only [lll|). 

Since light particles could yield gamma rays at energies 
that have been already probed experimentally, we mainly 
focus on their indirect detection signature to determine 
whether or not it is reasonable to consider them. We 
find that the gamma ray fluxes associated with particles 
lighter than 0(100 MeV) are in conflict with observations 
unless the v^-dependent term in the annihilation cross 
section i) is much larger than the S-wave term (devel- 
oped at the first order) and ii) satisfies the relic density 
requirement. Radio fluxes also validate this conclusion. 



Acceptable values of the cross sections 

Relic density calculations provide a strong constraint 
on any DM candidate. When DM particles are able to 
annihilate {i.e. when their non-relativistic transition oc- 
curs before their thermal decoupling), one obtains a sim- 
ple relationship, independent of mdm, between the DM 
cosmological parameter ^ll^h^ and the total annihilation 
cross section 12]. By requiring that fl*j!^h^ matches the 
observed value (^ 0.1, see [l^), one gets the following 
approximate annihilation cross section 
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with XF=mdm/TF ^ 17.2+ ln{g/ ln(md,„/ GeV) + 
\iiy/xF e[12-19] for particles in the MeV-O(GeV) range 
{g and g^^ being the number of internal and relativistic de- 
grees of freedom respectively). We write the cross section 
as avr ~ a -|- b{v/c)'^ where a and b are some constants 
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related to the S and P-wave terms (here v is the DM ve- 
locity and Vr the relative velocity). We shall assume first 
that b (and take c = 1). 

Eq.QJ sets the maximum value of {avr)ann one can 
use to compute gamma ray and radio fluxes from DM 
residual annihilations. Larger cross sections are possi- 
ble if Hdm 7^ ndm (but rcsidual annihilations are un- 
likely or even impossible since no, or few, anti-particles 
would be left after the DM freeze-out, tfo) or if the co- 
annihilation mechanism (which involves DM and X par- 
ticles) is at work In this case, however, the relation- 

^liiP ^ 1 + In {^y^' ^) indicates 
that mx should be close to rridm, which is actually ex- 
cluded if nidm < 0{ GeV) and if X is a charged particle. 
Thus, unless there exists a neutral (long-lived) particle 
X, the maximum annihilation cross section (times rela- 
tive velocity) into ordinary particles {e.g. 7, e~) that is 
legitimate to consider is about 10~^^ cm^^ s~^. 

The least massive annihilating WIMP 

Since eq. (Q) is almost independent of the DM mass (the 
only dependence in nidm being the logarithm which is 
"hidden" inxp), even light candidates are expected to be 
allowed. However, when dealing with the range mdm ^ 
O(GcV), having {avr)ann ~ lO^^^cm'^s"^ rather favors 
particles with an annihilation cross section independent 
of the DM mass or, for which the annihilations rely on the 
exchange of light neutral particles. Taken at face value, 
the first point strongly suggests that light DM should be 
made of scalars (which is in a gree ment with the second) 
and with non chiral couplings |1!tI | . One now has to check 
that their damping scale is not too large. 

In contrast with the free-streaming length [Ifs] of 
non-annihilating DM particles which depends only on 
TTidm, the scale Ifs of interacting particles in the 
[MeV,0(GeV)] mass ran ge d epends on both the DM 
mass and interaction rate [10|. By imposing If^ ^ 100 
kpc (the scale of the smallest galaxies) and assuming 
that interactions are weak enough to imply a decoupling 
in the radiation dominated era, one finds that particles 
with mdm ^ 0( MeV) (a^ec/lO"'') induce a cut-off in 
the matter power spectrum at ~ 100 kpc {adec being 
the scale-factor at the DM thermal decoupling). Nev- 
ertheless, in a realistic particle physics model, the ther- 
mal decoupling (based on the estimate of the DM-e elas- 
tic scattering cross section) is seen to occur around 1 
MeV {adec ~ 10-1°), so the mass range above 1 MeV 
should not modify the matter power spectrum at cosmo- 
logical scales. In fact, even if the DM thermal decou- 
pling was in the matter dominated era, only masses of 
mdm < 15 MeV (ildmh^ /O.l)-^ would affect the IO^Mq 
scale. The limit on mdm would actually get even smaller 
if DM was thermally decoupling after the non-linear col- 



lapse 01 (c./. self-interac ting DM for instance [l^ . 
which now appears unlikely [13 )j as the primordial fluc- 
tuations "disappear" to form objects. Collisional damp- 
ing due to J/-DM interactions may also contribute but 
is expected to yield a damping mass < IO^Mq so light 
annihilating candidates certainly deserve to be studied. 

We note that low DM masses are not constrained by 
direct detection experiments which are only sensitive to 
masses greater than ~ 7 GeV (except for e.g. Cresst 
[3, MACHe3, ROSEBUD, Tokyo ^ which are or will 
be able to go as low as ^ 1 GeV). In fact, exploring 
the low DM mass region should be a problem for cryo- 
genic detectors since the detection mechanism they cur- 
rently use (based on nucleus recoil) does not allow for the 
detection of particles much lighter than ^1 GeV with- 
out a significant effort (even by using the lightest pos- 
sible nucleus). Light scalars could escape searches in 
e+e^ colliders even when their production is based on 
the exchange of massive fermions (with a mass ^ 100 
GeV) as their cross section for anomalous single pho- 
ton events is still below (albeit very close) the sensitivity 
of past experiments 0| . Moreover, if they are able to 
annihilate into photons, the most interesting signature 
e+e^ e^e^l^l^ {Ip denoting the DM particles i.e. 
some missing energy) at the o!^ order should also be 
invisible because the electrons should remain mainly in 
the beam pipe, as one can infer from a kinematic anal- 
ysis. Nevertheless, even if there exists a deviation large 
enough to be detected, it is unlikely that PETRA or LEP 
experiments, for instance, got enough sensitivity. Hence, 
particle physics experiments still allow for the range we 
consider. Note that in any case DM should not have a 
coupling to the Z boson (otherwise it would have been 
detected in accelerator experiments) but since bino parti- 
cles or right-handed neutrinos, for instance, do not have 
this coupling either, this assumption seems reasonable. 

Indirect detection during the recombination epoch 

Let us check if DM annihilations at the recombina- 
tion epoch yield enough redshifted photons at an en- 
ergy g [keV, 10 MeV] to be detected nowadays. As- 
suming a ^ 5, one gets the following gamma-ray num- 
ber density w 2(cru,.)a„„ (j^dm)^ ^rec (with n^l^ 
and ircc the DM number density and Hubble time at 
the recombination epoch). This yields the present day 

flux $ - (^XiJ") (^«'->26 -m^^^ cm-2 sr^^ using 
$ = c 7i^/47r, {avr)26 — (o'Wr)ann/10-^^ cm-^ s^-^ and 
wmcV — "i(im/(MeV/c^). P-wave annihilation cross 
sections would decrease this flux "io"] but the latter 
is already overestimated since O(MeV) photons can 
lose some energy by scatterin g o n remnant free elec- 
trons and by reionizing atoms [2l|. Comparing $ with 

the much higher observed fluxes, namely ^[,fts'^*'''^°^(> 
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Emin) - 3 10-3(i;^i,/MeV)-i V(cm2ssr) and 
4(1.5) lO^VCcm^ssr) for the range [30, 100] MeV (and 
> 100 MeV respectively) we find that such residual 
annihilations are not ruled out by recent data. 
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TABLE I: Angular function F{9) and central 7-ray flux <!>(< 
1.5°) for difi'erent galactic DM profiles, Rsoi = 8.5 kpc and po 



chosen so that p{Rs, 



0.3 GeV/c^ cm-' [23 



Haloes are promising targets for DM annihilation 
radiation searches due to their high central den- 
sity and small cosmological distance '23^. Here, we 
assume a DM distribution of the form: p{r) — 

PO (r/r,)-^ [f + (r/r,)"]~^^"^^/" where (a,/?, 7) are 
given m Table n for the Navarro-Frenk- White (NFW), 
Kravstov et al. (KRA), modified isothermal (ISO) 
and Binney&Evans (BE) models. The number of pho- 
tons produced per unit of time is given by: n^(r) — 

no {p{r)/poY with no = 2 {(JVr)ann (Po/™<im)^- 

• The radiative flux from DM annihilations inside the 
Milky Way being distributed over a large fraction of the 
sky, we can compute the flux reaching the Earth from a 
circular cone with an opening angle 2 9 centered on the 
galactic centre: $(<6') = nor^F{0) = 6.2 IQ-^'cm-^s"! 
{avr}26rn^ly (po/(GeVcm-3))2 r,F (9)/ [10 kpc), where 
F{9) is a dimentionless function which depends on the 
DM profile. Numerical values of the flux are listed 
in Table □ They have to be compared with EGRET 
data from the galactic centre H^: $oL(< l-^"' ^7 > 
E^i^) = 2.2 10-6 cm-2 s-i {E^iJ{lOOMeV))-°-^ within 
[30 MeV, 1 GeV] (and even less within [1 GeV, 10 
GeV]) which is in agreement with the COMPTEL flux 
*oL(< 1-5°, > MeV) - lO-^cm-^s-i (or - lO^^ 
within 5°). 
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TABLE II: Expected fluxes from the Coma (C) and Virgo 
(V) cluster for different DM profiles |^. For the /3-profile of 
Virgo, only the flux within 1 Mpc is given, h — 0.7. 



• The total gamma ray flux from the DM halo of a 
galaxy cluster located at distance D is well approximated 
by $ci = JdVh-y{f)/{4TrD^). We list it for different halo 
profiles of two nearby clusters (Table (23|. These 
values have to be compared with <&o5s(> 100 MeV) < 
4 10^^ cm^^ s^^ at 2 sigma for these two clusters |2^ . 

By comparing the numbers in Table with 
we conclude that the annihilation cross section into 27 
(times relative velocity) of particles lighter than O(IOO) 
MeV should be much below (t7i'r)ann ~ 10-^^ cm"^ 3-^. 
This is clearly in conflict with the relic abundance re- 
quirement unless one assumes that i) avr is dominated by 
bv^ at the freeze-out epoch and ii) the a-term is smaller or 
equal to lO-^i - 10-^8 cm^ s-^ for m^m € [1, 100] MeV 
respectively (i.e. < 10^^ — 10^^ times the b-term at 
the freeze-out epoch). On the other hand, heavy parti- 
cles annihilating into photons and with nidm ^ 0{ GeV) 
appear quite compatible with observations even when 
aVr oc const (regardless of the profile used) . 

If the DM residual annihilations mainly proceed 
into leptons, then electrons could be detectable via 
synchrotron emission in a magnetic field B |27l |. 
For the galactic centre only, the flux Tu{< 0.1°) « 
5861.04 Jy(i//GHz)-0-5TOG^v {avr)26 (B/mG)-"-^ is 
compatible with the observed radio flux of SgrA (of the 
order of 360 Jy at = 330 MHz) provided nidm ^ 10 
GeV (assuming that all the annihilations proceed into 
e+e", a NFW proflle and B ^ 1 niG). However, masses 



< 10 GeV are OK if 



Majorana 



fermions that would annihilate into e'^e^ through the 
exchange of scalar particles may actually have this 
property (if they are chirally coupled and for some 
value of mdm) but they would fail to give the correct 
relic density if mdm ^ O(GeV). In contrast, scalars 
coupled to fermions have the correct relic density but 
not the correct radio flux. But these particles satisfy 
both criteria if they exchange a light gauge boson. For 
Coma (Virgo being swamped with emission from M87) 
only mdm 10 GeV produces radio emission at observed 
frequencies but the latter is (depending on B S [0.1, 10] p 
G) of the same order of magnitude as the observed flux 
for mdm ^ 10 GeV only (it is less for higher masses). 



Discussion and conclusion 

We now discuss three specific mass interval which have 
not been studied significantly before. They represent a 
"hole" in the investigations of the DM parameter space 
and deserve to be studied at least to probe whether or 
not there exist alternatives to heavy DM particles. 

nT'dm G [lMeV,m^[: Since the annihilation cross sec- 
tion of a pair of scalars into e'^e^ (via a fermion ex- 
change F and non-chiral couplings) is expected to be 
free of mdm (and potentially of the order of eq.JQ), see 
£-9- 113) J one can discuss the case of light scalar DM 
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candidates. However, to be compatible with the galac- 
tic centre COMPTEL/EGRET data, their annihilation 
cross section should be dominated by a term in a.ttfo- 
This is actually in conflict with the F exchanges (needed 
to get the correct relic density) but one can postulate an 
asymmetry between the DM and anti-DM number densi- 
ties or assume that annihilations mainly proceed through 
the exchange of either a light neutral fermionic WIMP or 
gauge boson. In any case, a more careful study is needed 
to ensure that the particles introduced satisfy all exper- 
imental constraints 0. Such light DM particles could 
perhaps be detected or excluded by using their interac- 
tions with ordinary matter or searching for the particles 
supposed to be exchanged in the annihilation process. 

fTT-dm G ['77.p,m7-[: They may be compatible with ob- 
served fluxes and relic density even for a ^ b but the 
production of D +^ He tends however to favor an 
annihilation cross section dominated by the term in w^. 

mdm G [m-r, 0(10 GeV)] : They seem in agreement 
with observations even when a ~ 10~^^cm'^s~^ but, if 
they mainly annihilate into e'^e~, radio fluxes tend to 
favor avr oc bv'^\tj,^. If they mainly annihilate into a 
pair T — f, one expects a soft gamma ray emission (no- 
tably from TT decays) plus an excess of positrons (due to 
f decays). Their energy is expected to be — {Ef — 

{Efsr)){1 + ^l-||cos0)/2 e [0,mdm - {Efsr)], 

where (Epsr) = ~ - m,/md^f+^ is 

the energy lost by the taus under the form of a sin- 
gle photon emission, i.e. via the Final State Radiation 
mechanism (FSR). As an illustration, with (3 ^ 0.09 
and mdm ^ 10 GeV, one finds (Efsr) ~ 1.5 GeV in- 
dicating that the emitted positrons could have an en- 
ergy of i?e+ — 8.5 GeV which seems of the order of the 
lower bound of the energy range indicated by HEAT 
experiment [2^. Of course, m^m ^ 0(10 GeV) will fail 
to explain the e+ excess above ~ 10 GeV and more care- 
ful studies are required to check whether particles with 
rrir < mdm < 0(10 GeV) can indeed escape present DM 
searches (which could perhaps be the case if they have 
stronger interactions with e~ than with nuclei, giving 
them a potential signature). 
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